Fractional Excitations in Cold Atomic Gases 
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We study the behavior of excitations in the tilted one-dimensional Bose-Hubbard model. In 
the phase with broken symmetry, fundamental excitations are domain-walls which show fractional 
statistics. Using perturbation theory, we derive an analytic model for the time evolution of these 
fractional excitations, and demonstrate the existence of a repulsively bound state above a critical 
center of mass momentum. The validity of the perturbative analysis is confirmed by the use of t- 
DMRG simulations. These findings open the path for experimental detection of fractional particles 
in cold atomic gases. 



Excitations carrying fractional quantum numbers are 
one of the most intriguing features of strongly interact- 
ing many-body systems. Arguably the most celebrated 
of those is the charge e/3 Laughlin quasiparticle respon- 
sible for the fractional quantum hall effect. Much effort 
is devoted to finding novel phases with even more ex- 
otic excitations such as non-abelian anyons and fractional 
statistics in three-dimensional systems. In the quest for 
the experimental realization of such strongly correlated 
phases, cold atomic gases with their clean and control- 
lable environment are a promising candidate and test- 
ing ground pQ. A first step in this direction marks the 
recently observed phase transition in a one-dimensional 
cold atomic system [2], where the fundamental excita- 
tions in the broken-symmetry phase are domain walls 
carrying an effective fractional charge. In this letter, we 
investigate the creation and detection of these fractional 
excitations in the one-dimensional tilted Bose-Hubbard 
model. 

A variety of theoretical proposals in cold atomic gases 
focus on the realization of strongly correlated phases with 
fractional excitations, such as models supporting spin liq- 
uid phases [3 , as well as Kitaev's toric code with abelian 
anyonic excitations [IHS], and systems in large effective 
magnetic fields [3 [8] exhibiting fractional quantum hall 
states [9HT5]. These proposals are based on standard 
experimental tools nowadays available in the context of 
quantum simulation, see Ref. [16 for a review. Espe- 
cially, the latest development of experiments [TTJ [18] with 
single-site readout and addressability in optical lattices 
has opened the path for the observation of a novel quan- 
tum phase transition in tilted optical lattices [2j [T9U2T] . 
The transition in this one-dimensional system takes place 
from a phase with one atom on each lattice site, to a 
broken-symmetry phase (BSP) with an alternating mean 
occupation, see Fig. [I] 

Here, we study the behavior of excitations in the 
tilted one-dimensional Bose-Hubbard model (TBH) in 
the phase with broken symmetry, with the latter be- 
ing two-fold degenerate. Fundamental excitations are 
domain walls between the degenerate broken-symmetry 
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FIG. 1. Phase diagram for the tilted Bose-Hubbard model 
in the regime A,w <C U, E > and an averaged density of 
one particle per lattice site: Mott insulator (MI) ground state 
for A/w < 7c ~ 1.85 [19] , and the broken symmetry phase 
(BSP) for A/w > 7c . 



phases and exhibit an effective fractional charge. In turn, 
the simplest excitation corresponds to moving a single 
atom from one site to its neighboring site. Such an exci- 
tation corresponds to the creation of two closely bound 
domain walls. The important question is then, whether 
these experimentally accessible excitations will decay into 
the fundamental domain-wall excitations and how to de- 
tect these fractional excitations. We derive an analytic 
expression for the time evolution of those excitations, 
and show the existence of a repulsively bound state of 
fractional excitations above a critical center of mass mo- 
mentum Q c = 27r/3. Further we provide experimental 
signatures for measuring fractional excitations in a setup 
of cold atoms and give direct numerical calculations for 
the time evolution for finite size sample. 

We start with the one-dimensional tilted Bose- 
Hubbard model as realized in the experimental setup [17] , 
which takes the form 

H TBH = - W ^ ( a i a *+l + a i a \+l) 

+ 77 / ] ni(ni - l) - Eyjm 

i i 

with the on-site interaction J7, the hopping rate w, and 
the lattice tilt E per site. In addition, a\ (o^) is the cre- 
ation (annihilation) operator for a particle on site z, and 
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FIG. 2. Mapping from the tilted Bose-Hubbard model 
(TBH) onto the effective hard-core boson model (HCB): two 
neighboring lattice sites occupied by a 0, 2-configuration map 
onto a single unoccupied site, whereas a singly occupied site 
maps onto an occupied site. A second-order transition with 
2w 2 1 A allows for an effective hopping of the fractional exci- 
tations. 

Ui = a\di is the number operator on site i. We focus on 
the regime with an averaged density of one particle per 
lattice site, and assume a positive tilt E > 0, i.e. the lat- 
tice is tilted to the right. The system is in a meta-stable 
state, where the condition that the on-site interaction U 
and the lattice tilt E are much larger than the tunneling 
energy w, prevents the relaxation into a state with more 
than two-particles accumulating in a single site. Then, 
the relevant energies are the energy difference A = E — U 
and the hopping rate w. 

For large negative A, the ground state is a Mott insu- 
lator (MI) with one particle per lattice site. Increasing 
|A| ~ w, particles can tunnel to their right neighboring 
site, as long as the particle on that site has not tunneled. 
Eventually, the system undergoes a phase transition into 
a ground state with broken translational symmetry with 
two atoms on each second lattice site, see Fig. [I] 

The condition E,U ^> w, A suppresses processes with 
particles hopping to the left as well as the accumulation 
of more than two particles in a single lattice site. These 
constraints are most conveniently incorporated by a map- 
ping of the Bose-Hubbard model to a spin model: The 
spin degree of freedom resides between two lattice sites 
i and i + 1 and the spin-up state corresponds to a par- 
ticle at lattice site z, while the spin-down state accounts 
for the situation where the particle has tunneled to the 
site i + 1. Then, the Bose-Hubbard model maps to a 
one-dimensional Ising model 

if Ising = - 2V2wJ2si + A^>* 




with S\ (S l x ) is the spin operator along the z (x) axis. 
In addition, the last term accounts for the constraint 
that a particle can only tunnel from site i to i + 1 if 
there is already a particle at site i + 1 and formally 
requires taking the limit J —> oo. Under this map- 
ping, the Mott-insulating state corresponds to a param- 
agnetic phase with all spin- up, while the ordered broken- 
symmetry phase corresponds to an Ising antiferromag- 
netic ground state. 

In the following, we focus on the quantum phase for 
A ^> w with two particles on every second lattice site. 



The experimentally accessible excitations are achieved by 
moving a particle from one lattice site to its neighboring 
site to the left, which corresponds to a single spin- flip. 
The energy of this excitation is given by A. This single- 
particle excitation can now be decomposed into two frac- 
tional excitations via a second order-process as indicated 
in Fig. [2] In the classical regime with w = 0, these delo- 
calized fractional excitations have the same energy A as 
the single-spin excitations, while adding quantum fluc- 
tuations within second-order perturbation theory yields 
a nearest-neighbor interaction V = 2k; 2 /A, as well as 
an effective hopping for the fractional excitations with 
amplitude z = 2w 2 / A. Note that fractional excitations 
are restricted to even or odd lattice sites, depending on 
the local order, and hence hopping takes place from site 
i to site i ± 2. This restriction can be uniquely fulfilled 
by introducing a new lattice: each singly occupied site 
and the combination of a doubly occupied site next to 
an empty site count as a new lattice site (see Fig. [2|. 
Consequently, the number of lattice sites between two 
excitations is halved, and the dynamics of fractional exci- 
tations is governed by the hard-core boson model (HCB) 
with nearest-neighbor hopping and interaction z and V, 
respectively, described by the Hamiltonian 

H HCB = ~z^2 ( b ] b J+i Jrb 3 b )+i) + V ^2 n 3 n 3+i' ( 3 ) 

3 3 

Here &t (bj) is the creation (annihilation) operator for 
a fractional excitation on site j, and nj = b^bj, is the 
number operator with (nj) G [0, 1]. 

A similar Hamiltonian has previously been studied in 
the context of repulsively bound pairs [22j [23] . Here, 
a bound state of domain walls means that the experi- 
mental excitation cannot decay into its fractional parts 
due to energy conservation, while the absence of a bound 
state corresponds to a situation where this experimental 
excitation decays into delocalized fractional excitations 
and allows for their experimental observation. As in the 
present situation, z = V is the only energy scale of the 
model, we obtain a universal bound-state structure for 
the fractional excitations, which we analyze in the fol- 
lowing. 

We write the wave-function for two fractional excita- 
tions in the effective model as ip(i,j), where i and j de- 
note the positions of the fractional excitations in the ef- 
fective lattice. The discrete translational invariance of 
the system provides conservation of the center of mass 
quasi-momentum, which allows for an expansion into 
eigenfunctions of H HCB for fixed center of mass quasi- 
momentum Q. Therefore, the two-particle wave-function 
can be factorized as ip(X : x) = e~ z ® x \jj(x) with the cen- 
ter of mass X = (i + j)/2 and relative coordinates x = 
i — j. We find two different regimes for the two-particle 
states: In the first regime for \Q\ < Q c = 27r/3, the two- 
particle eigenfunctions are given by scattering states ifj q 
alone with energy E q = — 2zq cosg, where q is the rela- 
tive momentum and zq = 2z cos Q/2 denotes the hopping 
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rate in the center of mass frame. Its wave- function re- 
duces to plane waves, ip q (x) = (1 — 5 Xj q) cos(q\x\ + </>Q, g ), 
with a scattering phase shift 
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The Kronecker-delta factor 1 — accounts for the hard- 
core constraint and enforces ip q = at x = 0. In the 
second regime for \Q\ > Q c , an additional bound state 
iJj b emerges. The repulsive interaction yields an energy 
E B = z(l — 4 cos 2 Q/2) lying above the scattering contin- 
uum. Its two-particle wave-function shows an exponen- 
tial decay with relative distance x, and can be written 
as 



i/j b (x) = (1 - ^,o) 



-4 cos 2 Q/2 
4 cos 2 Q/2 



-2 cos 



1*1 



Note that the alternating amplitude of the wave-function 
is a typical feature of a repulsively bound state. The 
general wave-function of two fractional excitations with 
fixed center of mass momentum Q can therefore be de- 
composed as 



-iE B t/h 



I 



Cqt/Jqe 



-iE q t/h 



(4) 



with C B and C q denoting the overlap of the initial wave- 
function with the bound and scattering eigenfunctions of 
H nC B- The experimentally accessible initial states are 
achieved by moving one particle to the left, which in 
the effective model corresponds to two adjacent occupied 
sites ^spe(^) = S Xj i. Then the overlap is given by C B = 
6(\Q\ - Q c )v/l-4cos 2 Q/2 and C q = V2cos(q + Q>g ), 
with Q(x) being the Heaviside step function. The inte- 
gral over the relative quasi-momentum q in eq. Q can be 
carried out analytically, giving rise to a formal solution 
in terms of an infinite sum of Bessel-functions of the first 
kind 



^C q il> q (x)e- iE *«* 



Y,c n {x)e^ n J n {2tz Q /K) (5) 



with coefficients c n (x) defined via a discrete Fourier 
transform C q ip q (x) = ^Z n c n (x)e zqn . 

The time evolution of the wave-function is shown in 
Fig. [3j The superposition of the scattering states leads 
to a ballistic expansion of the fractional excitations with 
a velocity determined by the hopping energy zq, i.e., 
the two-particle wave-function ^ SPE is strongly localized 
around a linearly growing relative distance x ~ tzQ/h, 
with some additional interference fringes appearing at 
smaller relative distances x, but propagating at the same 
velocity. However, the finite overlap with the bound state 
ip B for Q > Q c creates an additional stationary peak 
at x ~ 1 (see Fig. |3jb)). With the scattering states 
moving away from each other, measurement of the wave- 
function amplitude at x = 1 at times t ^> 1/zq allows 
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FIG. 3. Time evolution of the relative wave-function ip(x): 
(a) For \Q\ < Q c , the wave-function is strongly localized at 
around a relative distance x ~ tzQ/h. (b) For \Q\ > Q c , 
the bound state adds an additional exponentially decaying 
contribution, with its maximum at x — 1. 



one to single out the bound-state contribution. For- 
mally, this can be cast in terms of a correlation function 
C(t) = |^ SPE (X, x = 1)| 2 on the effective lattice. In the 
microscopic lattice, C(t) corresponds to (P^P^i), with 
Pi = rii(ni — 2) being the projection operator on singly 
occupied sites. The time evolution of C(t) for different 
center of mass quasi-momenta is shown in Fig. [4j For 
Q < Qc the correlation function decays to zero with a 
characteristic behavior ~ (tzQ/h)~ 3 , and exhibits char- 
acteristic oscillations due to interference between the dif- 
ferent scattering states. In addition, the decay exhibits 
an intermediate regime with a characteristic behavior 
~ (tzQ/h) -1 . The time scale for the crossover between 
the long-time behavior and this intermediate regime di- 
verges approaching the critical value Q c . Consequently, 
the decay at Q = Q c is given by a critical behavior 
~ (tzQ/K) -1 for the correlation function, which follows 
from the analytical expression for the wave-function 

^peW =e-^M [J lx] (2zt) +iJ\ xhl (2zt)] . 

Finally, the presence of a bound-state is characterized by 
a saturation of the correlation function 

C(t)^0(\Q\-Q c ) [l-4cos 2 Q/2] 2 

for tzQ/H ^> 1. The bound-state contribution grows 
steadily towards the edge of the Brillouin zone, and even- 
tually, V^spe becomes an eigenstate of the Hamiltonian at 
\Q\ = 7r with a constant correlation function C(t) = 1. 
In addition, the interference between the bound state 
and the scattering states gives rise to a characteristic 
beating of the correlation function with the frequencies 
cj± = E B ±E sc (q = 0) = z(l -4cos 2 Q/2±2cosQ/2). 
This gives experimental access to the energy gap between 
bound and scattering states, and allows for the measure- 
ment of the bound-state energy. 

To confirm our perturbative results and to provide 
further insight into an experimental realization, we pro- 
vide time-adaptive density matrix renormalization group 
(t-DMRG) simulations [24l [25] of a single-particle ex- 
citation in the original tilted Bose-Hubbard model in 
eq. ([T]). Here, we use a realistically large lattice of L = 30 
sites, i.e. I = 16 sites in the effective lattice. Time- 
evolution calculations are performed using second- and 
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FIG. 4. Time evolution of the correlation function C(t) for 
different center of mass quasi-momenta Q: Below the critical 
momentum \Q\ < Q c , the correlation function decays to zero, 
while for \Q\ > Q c , the overlap with the emerging repulsively 
bound state gives a finite probability for fractional excitations 
to stay at finite relative distance x, resulting in a finite value 
of the correlation function for tzQ/h ^> 1. 



fourth-order Trotter decompositions. A comparison be- 
tween the analytic correlation function and the t-DMRG 
result is shown in Fig. [5] 

On the one hand, t-DMRG results show good agree- 
ment with the perturbative model. We find that the 
timescale of the correlation-function decay and the satu- 
ration values agree well with the analytic results. 
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FIG. 5. t-DMRG results (solid lines) in a finite lattice with 
L — 30 microscopic sites (i.e. / = 16 sites in the effective 
lattice) in comparison to the analytic results (dashed lines) 
in an infinite lattice for various center of mass momenta Q. 
Finite size effects start come into play at times t ~ Ah/ z. 

On the other hand, we observe deviations due to fi- 
nite system size and finite values U, E, and w, which we 
discuss in the following. First, the simulations are per- 
formed in a finite size system. This leads to revivals of the 
correlation function C(t) due to scattering of fractional 
excitations at the system boundaries. In a microscopic 
lattice of L = 30 sites, we find deviations due to finite 
size at times t ~ 4h/z. A hard-core model neglecting 
nearest-neighbor interaction, which allows for an analytic 
solution even in a finite system, agrees well with the t- 
DMRG result and gives an estimate for revival times. 
From this, we can derive a lower boundary of L > 26 



to observe a signature for a bound state in the correla- 
tion function. Second, the simulated correlation function 
shows an additional high-frequency oscillation. This can 
be explained as follows: Since the initial state is not an 
energy eigenstate of the Hamiltonian for a finite ratio 
w/A, this leads to interference between different energy 
contributions, and ultimately in oscillations on the order 
of the excitation spectrum, i.e. A. We find that for a 
ratio of A/w > 30 the suppression of these oscillations is 
strong enough to see a clear correlation-function signal 
(see Fig. |5|. 

Finally, the validity of the correlation function is based 
on the stability of the broken- symmetry phase. However, 
second-order processes in w/U and w/E allow the sys- 
tem to reach states with more than two particles on a 
single lattice site (see [H]). Our simulations show that 
the probability for having three particles on a single lat- 
tice site at times t = 4H/z of the order 10 -6 . It follows 
that this effect can be safely neglected for an on-site in- 
teraction U > 8A on experimentally relevant timescales. 
With typical lifetimes of atoms in an optical lattice be- 
ing on the order of seconds, an effective hopping rate 
z/h ~ 1 Hz should be sufficient for observing the satura- 
tion of the correlation function C(t). Further assuming 
a microscopic hopping rate of w = 15 Hz, suppression of 
high-frequency oscillations yields a single-particle excita- 
tion energy A/h = 450 Hz. Finally, an on-site interaction 
U/h = 3.6 kHz results in a stable broken-symmetry phase 
for the duration of the experiment. 

An important aspect is that throughout this paper cal- 
culations used an initial state of delocalized excitations 
with finite center of mass momentum Q in order to max- 
imize the correlation function C(t). This behavior may 
seem counterintuitive, as one would expect a localized ex- 
citation to yield stronger correlations. However, localiza- 
tion in configuration space gives rise to a flat momentum 
distribution; with C(t > H/z) -» for \Q\ < Q c , this 
averaging over center of mass momenta then results in a 
decrease of the correlation- function by a factor ~ 0.17. 
Still, the finite-Q initial state can be composed of local- 
ized excitations via superposition, i.e 

IV>Q> = E e^ j al_ iaj \BSP) (6) 

j'G2N 

with |BSP) denoting the broken-symmetry phase with 
two particles on every even lattice site. The factor of two 
in the phase then accounts for the two-site hopping of 
fractional excitations in the microscopic model. 
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